Density expansion for the mobility in a quantum Lorentz model 
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We consider the mobility of electrons in an environment of static hard-sphere 
^ ■ scatterers, which provides a realistic description of electrons in Helium gas. 

A systematic expansion in the scatterer density is carried to second order 
relative to the Boltzmann result, and the analytic contribution at this order 
is derived, together with the known logarithmic term in the density expansion. 
It is shown that existing experimental data are consistent with the existence of 
the logarithmic term in the density expansion, but more precise experiments 
are needed in order to unambiguously detect it. We show that our calculations 
provide the necessary theoretical information for such an experiment, and give 
a detailed discussion of a suitable parameter range. 

PACS numbers: 51.10+y, 05.60+w 
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I. INTRODUCTION 



It has been known for almost thirty years that transport coefficients, in contrast to 
thermodynamic quantities, do not possess a virial expansion.0 This insight came originally 
as a substantial surprise,^ and the physics behind it turned out to be very fundamental 
in nature. Namely, it is related to the, previously unexpected, existence of long-range 
dynamical correlations in equilibrium fluids. As a specific example, let us consider the 
diffusion of a tagged particle in a fluid. The classic method for setting up a virial or 
density expansion for the diffusivity was based on a generalized Boltzmann equation.! In this 
method one constructs a virial expansion for the collision operator by taking into account 
collisions of the tagged particle with a successively increasing number of scatterers which 
is, however, finite at every order in the expansion. This is analogous to the Ursell-Mayer 
cluster expansion for thermodynamic guantitiesE Alternatively, one can use the Green-Kubo 
or time correlation function method.0 Both methods yield identical formal results for the 
density expansion of the diffusivity D (or any other transport coefficient), 

D/D B = 1 + Din + D 2 n 2 + 0(n 3 ) , (1.1) 

where D B is the Boltzmann result for the diffusivity, n is the dimensionless density of the 
fluid, and Di, D 2 , etc. are the 'virial coefficients'. The surprise consisted in the realization 
that the coefficients in this expansion are infinite past a certain order that depends on the 
dimensionality of the system.! It was soon realized that these infinities signalize the presence 
of a logarithmic density dependence in the expansion given by Eq. ( |1 . 1| ) B For instance, in a 
three-dimensional (3-d) system, the coefficient D 2 is n-dependent, and goes for small n like 
D 2 ~ Inn + const, so that Eq. (|1.1|) must be replaced by, 

D/D B = l + D in + D 2log n 2 \nn + D 2 n 2 + o(n 2 ) , (1.2) 

where D 2 i og and D 2 are numbers, and o(n 2 ) denotes terms that vanish faster than n 2 for 
n —>■ 0. The functional form of the higher coefficients is not known. The physical reason 
for the breakdown of the density expansion, and the appearance of the logarithmic terms, is 
a collective effect: Ring collisions, in which the scattered particle collides with a scatterer, 
and then visits a number of other scatterers before returning to the original one, turn out 
to lead to infinite phase space volumes if a finite number of scatterers is considered, and 
hence the distances between scatterers, and the times traveled between them, are allowed 
to become arbitrarily large. These collision events are, however, unphysical, because the 
inevitable presence of scatterers that do not belong to the particular cluster under consid- 
eration does not allow the effective distance traveled between successive scattering events 
to be substantially longer than a mean-free path. Properly taking into account this mean- 
free-path damping of the scattered particle's trajectory leads to the mentioned logarithmic 
term. The latter represents dynamic correlations in the fluid which extend over a distance 
of a few mean-free paths (or a time equal to a few mean- free times), which in the limit of 
small densities becomes much larger than the effective particle size, or the range of the in- 
termolecular force. In contrast, static correlations in classical equilibrium fluids extend only 
over the range of the intermolecular force, which is the reason why the cluster expansion 
works for thermodynamic quantities. 
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Due to its bearing on fundamental aspects of the Statistical Mechanics of both equi- 
librium and non-equilibrium systems, and its close connection to related effects, such as 
long-time tails,0 the breakdown of the density expansion for transport coefficients has gen- 
erated substantial interest over almost thirty years. The existence of the logarithmic term 
has been ascertained theoretically for a variety of classical systems, both real fluids and 
model fluids, and for a variety of transport coefficients.! It has also been seen in computer 
simulations of a 2-d Lorentz gas,i where the logarithm appears at first order rather than at 
second order due to the lower dimensionality. It has also been established that the same 
effect occurs for quantum mechanical scattered particles.i However, its experimental verifi- 
cation has proven to be extraordinarily difficult.il The reasons for this are manifold. First, 
the detection of any logarithmic term on an analytic background is very difficult due to the 
slow variation of the logarithm. Second, for 3-d systems the effect appears only at second 
order in the density expansion beyond the Boltzmann equation. Third, the coefficients in 
Eq. ( |1.2| ), or in the equivalent expressions for other transport coefficients, are not known 
for classical fluids with realistic interaction potentials. Even for a hard-core model fluid it 
has so far not been possible to calculate the coefficient D2 of the analytic term at second 
order. Since the detection of a logarithmic term on an unknown background is a hopeless 
task, this essentially precludes the use of 3-d classical fluids for a convincing observation. 
On the other hand, 2-d classical systems are very hard to realize. Furthermore, estimates of 
the values of the coefficients in classical systems show that for those transport coefficients 
that can be accurately measured, like e.g. the shear viscosity, the coefficient of the analytic 
second order term is much larger than the one of the logarithmic term. This is because the 
former contains essentially excluded volume effects, which dominate over the ring collision 
contributions. 

These difficulties raise the question whether quantum systems are possibly better suited 
for an experimental verification of the logarithmic term than classical ones. 2-d quantum 
systems, although easy to realize, are not suitable since transport in 2-d quantum systems 
is pathological due to localization effects. 0111 This leaves 3-d quantum systems. Here the 
effect also occurs only at second order, but otherwise the situation is much more favorable 
than in the case of 3-d classical systems, as we will see. A chief advantage is the fact that 
the hard-core Lorentz gas model, i.e. a tagged particle moving in an array of static hard- 
sphere scatterers,0 is a much better approximation for certain quantum systems than for 
any classical ones. 3-d quantum systems are therefore the most promising candidates for an 
experimental observation of the nonanalyticity. 

A particularly promising system consists of electrons injected into Helium gas of density 
n.0 The electron-Helium scattering process is well known, and its characteristics are conve- 
nient from a theoretical point of view. The scattering length, a s = 0.63A, is positive, and for 
thermal electrons the energy dependence of the scattering cross section is negligible. Since 
the electrons behave quantum mechanically, the thermal wavelength, A = (2ir 2 h 2 p/my/ 2 , 
with (3 = 1/kgT and m the electron mass, provides an additional length scale besides a s and 
the mean Helium atom separation n~ 1//3 . The leading parameter in the density expansion 
is naz X = A/4tt/,§ with / = l/Anna 2 the mean free path, and a s /X serves as an additional 
small parameter. The mass ratio m^ e /m pa 10 4 makes it a good approximation to treat the 
Helium atoms as static scatterers. Finally, the low density of the injected electrons allows 
one to neglect Coulomb interaction effects between the electrons. An experiment^ which 



3 



measures the mobility of the electrons (which is related to the diffusion coefficient by an 
Einstein relation) thus constitutes an almost ideal realization of a 3 — d quantum Lorentz 
model. 

The density expansion for the transport coefficients of the 3 — d quantum Lorentz model 



has been considered in Refs. fRlO. The experimentally relevant quantity is the mobility, /x, 



at finite temperature. The leading terms in the expansion for /x, analogous to Eq. (|1.2|) , are, 
fi/fia = 1 + ViX + V2iogX 2 ^x + H2X 2 + 0(xa s /X) + o(x 2 ) , (1.3a) 

with 

^ = -^2/6 , (1.3b) 



IHUog = ~ 4)/32 . (1.3c) 

Here /xb(T) = (4eZ/3)(27rm//3) _1//2 is the Boltzmann mobility, and x = A/7rZ and a s /X are 
small parameters. 

Adams et al.0 have used Eqs. ( |1.3| ) to analyze experimental data obtained from time- 
of-flight measurements for electrons in He and H 2 . Their main objective was to refute the 
popular misconception that /xi = /x 2 ; og = 0,0 which arose from an inappropriate application 
of localization ideas to the low-density regime. Ref. concluded that the existing experi- 
ments give very good agreement with the value of /xi given in Eq. (|1.3b| ). This success raised 
the question whether the same system could be used to observe the logarithmic term. In 
the absence of information about /x 2 this would involve measuring the conductivity over a 
gas density range that is sufficient to observe the logarithmic dependence directly. This is 
clearly hopeless. However, if /x 2 was known, then the logarithmic term would just provide a 
weakly density dependent correction to it, and a sufficiently accurate experiment at fixed gas 
density would be sufficient to probe the existence of the logarithmic term. This background 
provided the motivation for a calculation of /x 2 , the result of which has been reported in a 
recent short communication as," 

/x 2 = 0.236... . (1.3d) 
The purpose of the present paper is to provide the technical details of the calculation 



whose result was reported in Ref. [19]. We will also give a more detailed discussion of how 
our result can be used to design an experiment capable of observing the logarithmic term. 
The paper is organised as follows. In Section |J we first define the model, and set up 
a diagrammatic perturbation theory for the conductivity at zero temperature. We then 
perform the calculation to second order in the scatterer density, and convert the result into 



an expression for the experimentally relevant temperature dependent mobility. Section [II I 
contains a detailed discussion of the relevance of our results for a proposed experiment to 
detect the logarithmic term in the density expansion. Some technical details are relegated 
to two appendices. 
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II. DENSITY EXPANSION FOR THE MOBILITY 



A. The Model 

We consider a model of noninteracting electrons moving in three dimensional space, 
and scattering off static, uncorrelated, and randomly located impurities with density n. The 
electron-impurity interaction can be characterised by the scattering cross section a s = 47ra^, 
where a s is the scattering length. 

The Hamiltonian for this system is, 

^ = H(ek-^)4, ( x a k, ( x+ H ^(q)4,a«k-q,a , (2.1) 
k,<r k,q,(j 

where a k a , a k(J denote creation and annihilation operators for the electron with wavenumber 
k and spin a, fi is the chemical potential, ^(q) is the Fourier transform of the electron- 
impurity scattering potential, and e k = k 2 /2m is with m the electron mass. In this section 
we use units such that % — 1. 

We will use the standard Edwards diagram technique. Accordingly, we consider re- 
tarded (R) and advanced (A) zero-temperature Green's functions, 



£k, P M = 

and their impurity averaged counterparts, 



1 



u - H ± zO 



P (2.2a) 



= 4, P — \ t r,a ( \ > ( 2 - 2b ) 

where e F = fi(T = 0) is the Fermi energy, and E^'^u;) is the self energy, u measures the 
energy distance from the Fermi surface. Ultimately we will take the limit u —>■ to obtain 
the static mobility or conductivity. 

We will restrict our considerations to pure s-wave scattering, which is equivalent to 
the assumption of point-like scatterers. This approximation will substantially simplify all 
calculations. Its justification and limitations will be discussed in Sec. |TJ. For pure s-wave 
scattering the self-energy reads to linear order in the impurity density, 

£k' A (u;) = A ± i/2r + 0(n 2 ) , (2.3) 

with t = Im/kF the mean-free time. The real part of the self energy, A, is to this order 
independent of the wavenumber k. Even though it is formally infinite for pure s-wave 
scattering, it therefore strictly renormalizes the chemical potential or the zero of energy, 
and can be neglected. Although we will want to effectively expand the full Green's function 
to second order in the impurity density, we will find it convenient to use the self energy as 
given by Eq. (|2.3| ) and to include higher order contributions to the self energy explicitly at 
a later stage. We therefore define, for later reference, an auxiliary zero frequency Green's 
function, 

G? A = — . (2.4a) 
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We will also need the free electron Green's function, 

G (0)R,A = 1 (24b) 

B. Diagrammatic Expansion of the Conductivity 

It is convenient to calculate the conductivity, a, of degenerate electrons at T = 0, and 
then to convert to the experimentally relevant finite- T mobility by means of an Einstein 
relation and a Kubo-Greenwood formula. A general scheme for a diagrammatic calculation 
of the conductivity within the context of the model defined in Sec. [II A| above has been given 



in Ref. [T5|. The result was a density expansion for the conductivity of the form, 

a /a B = 1 + * x -L + O2log (J-)\ n (_L) + ff2 (_LV + (a s /l) + o (l/(M) 



(2.5) 



with <tb the Boltzmann conductivity, /cf the Fermi wavenumber, and a s kp considered small. 
There is no need to repeat the description of the general scheme here. Suffice it to say that 
the real part of the dynamical conductivity of noninteracting electrons at zero temperature 
can be expressed in terms of the Green's functions defined in Eq. (|2.2j ) by means of the 
Kubo-Greenwood formula, 

Rea{uj) = A_ Re ^ V (k) (g£ J ,(u)g&(u = o)-g^{u)G^(u = o)) «( P ) , (2.6) 

where v(k) — k • q/g, with q an arbitrary fixed vector, is the current vertex. The brackets 
(• • •) denote the averaging over the random positions of the impurities. We will perform this 
average by means of standard diagrammatic methods.0 We note that this is not the only 
possible way to calculate the conductivity. One could, for instance, use quantum kinetic 
theory instead. However, it is known from calculations of the coefficient o\ in Eq. (|2.5|) 
that the diagrammatic methodEl'B leads to much simpler calculations than quantum kinetic 
theory.0 Technically, our calculation is a systematic extension of the work of Refs. O and 



|16| . These authors had already identified all diagrams that contribute to 0\ and a 2 iog, and 
calculated both coefficients. For our purposes we will now have to calculate the previously 
considered diagrams to 0(n 2 ) instead of 0(n 2 In n), and to identify and evaluate all diagrams 
that contribute to 02, but neither to o\ nor to a^iog- 

Let us start by introducing symbols and abbreviations for the diagrammatic elements 
of our perturbation theory. The averaged Green's function, Eq. ( |2.4a| ), we denote by a 
directed full line, with an arrow pointing right and left for the retarded and advanced Green's 
function, respectively. We further denote the difference between the Green's function given 
by Eq. ( [2.4a|) and the free electron Green's function, Eq. ( |2.4b|) , by a directed line which 
carries a triangle. The impurity potential is denoted by a dashed line, and the impurity 
density by a cross; each cross with a dashed line running through it corresponds to a factor 
u = 1/2ttNft, with Np = kprn^TT 2 the free electron density of states per spin at the 
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Fermi level. The current vertex v (k) is represented by a triangle. All of these diagrammatic 
elements can be seen in Fig. [l]. We also define the quantities e = 2mep 1 and 7 = m/r. 

In terms of these quantities, the simplest diagrammatic contribution to Eq. ( |2.6| ) in 
the limit of zero frequency is the simple bubble shown in Fig. 1(a), and the corresponding 
diagram with the direction of the lower line reversed. The analytic expression for this 
contribution to the static conductivity, which we denote by o"(i a ), reads, 

„2 



O-(la) 



Re £ Kk) 



Performing the integral is simple, and yields, 

3 ' ~ n 2 



°"(la) — &B 



1 + 



+ 




(2.7a) 



(2.7b) 



Diagram 1(a) thus contributes 3/2 to the coefficient 02, a result which we record in Table |. 
For later reference we also give explicitly the Boltzmann result for the conductivity, 



e 2 e 



3/2 



37T 2 7 



e 2 n e r 
m 



(2-7c) 



with the free electron density n e = e 3//2 /37r 2 , and e the electron charge. Notice that ob is a 
function of e, this will be of importance later. We also note that the simple bubble is the 
only diagram for which we have to consider the Q R Q R contribution to Eq. (|2.6|) . All other 
diagrammatic contributions to this term turn out to be of higher than second order in the 
impurity density. 

Figure [I] also shows two diagrams that have to be considered due to our choice of the 
basic Green's function, Eq. (|2.4a| ). If we had calculated the self energy in the basic Green's 
function to second order in the impurity density these diagrams would not appear. Note 
the appearance of the 'triangulated' Green's function, which is necessary to avoid double 
counting. The evaluation of these diagrams is also straightforward, and the results are given 
in Table |. 

We now turn to the diagrams that were discussed previously in Refs. ^,0, and were 
calculated there to 0(n 2 lnn). They are shown again in Fig. |^. We find that all of these 
diagrams contribute to the analytic term at second order as well. The diagram shown in 
Fig. 2(a) is of special interest, since for pure s-wave scattering it is found to be ultraviolet 
divergent. We therefore introduce an ultraviolet cutoff Q ~ l/a s , and discuss this diagram 
in some detail. We stress that this divergence is due to an unphysical treatment of the short- 
range part of the electron-impurity interaction and often occurs when one uses an s-wave 
scattering approximation. It should not be confused with the physically more interesting 
logarithmic singularity in the density expansion of the transport coefficients which is due to 
long-range collective effects. 

The analytic expression corresponding to Fig. 2(a) is, 



, 7 



a (2a) = 2u 2 £ Kk)] 2 Re (G r G 

k,p,q 

Re / dqq 2 J ++ {q) 



R (~iR 
k p ^k-p- 



(~iR riR (~iA 
q L7 k-q ^k °k 



2vr 5 e 3/2 



(e-ry) \J ++ {q) ~ J + ~ 
-(e + z 7 ) 7 f- J ++ (q) 



) 

(<?)] 



(2.8) 
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To obtain the second equality we have repeatedly used the identity 



G k G t - — ( G k - G k) > (2-9) 



and defined two functions, 



J + -(q) = / dk — — ^— , (2.10a) 

J e — k 2 + 27 e — (k — q) z — ry 

and 

J ++ (q) = [ dk -A— \ . (2.10b) 

J e — k 2 + 27 e — (k — qj 2 + ij 

The same functions had been defined in Refs. 0,0, and all diagrams can be expressed 
in terms of them. For our present purposes we need a more accurate evaluation of these 
functions than the one that was given before. It turns out that both integrals can be done 
exactly in closed form, a task which we relegate to Appendix |A]. Using the result in Eq. 
a calculation sketched in Appendix [B] leads to, 



C(2a) = °~B 



-7T 



(2.11) 



The In Q-contribution stems from the real part of the self-energy contribution to the upper 
electron line in diagram 2(a). It thus constitutes just a shift (albeit an infinite one in the 
case of a pointlike potential) of the chemical potential, like the constant A in Eq. (|2.3|) which 
we neglected earlier. We therefore expect the InQ-term to disappear upon considering the 
experimentally relevant mobility instead of the conductivity. We will find this expectation 
to be borne out later, cf. Sec. [II C| . The remaining diagrams in Fig. |2| can be calculated 
along the same lines, and the results are given in Table |. 

The diagrams shown in Figs. 2(d)-(h) allow for generalizations which all contribute to o<i- 
These generalizations are obtained by replacing the 'ladder' and 'crossed-ladder' elements in 
these diagrams by the respective infinite resummations. The resulting diagrams are shown in 
Fig. |. Note that the ladder, or 'difluson', and crossed-ladder, or 'Cooperon', resummations 
in Fig. |H start with three rungs each to avoid double counting of the diagrams of Fig. |2j. The 
infinite resummations are again easily expressed in terms of integrals over the functions J ++ 
and J + ~ , and the results are listed in Table |. The reason why these diagrams with increasing 
numbers of impurity lines, and hence increasing numbers of factors 7 ~ n, all contribute to 
the same order lies in the fact that the diffusion pole contained in the ladder and crossed- 
ladder resummations leads, with increasing order, to increasingly singular infrared behavior 
of the integrand, which is cut off only by 7. The cancellation of these two effects leads to all 
of these diagrams being of the same order in 7 or n. By the same argument it follows that 
these are the only infinite resummations (given our definition of skeleton diagrams by means 
of Eq. ( p.4a| )) that contribute to the desired order. In particular, diagrams that contain 
more than one diffusion pole do not contribute. 

We now turn to other skeleton diagram contributions. All relevant diagrams with three 
impurity lines (as far as they were not included in Fig. |2|) are shown in Fig. |], and all those 



S 



with four impurity lines are shown in Fig. |5|. The evaluation of these diagrams offers no 
particular difficulties, except that some care has to be exercised since the vector nature of 
the current vertex leads to some nontrivial angular integrations. Of course, all diagrams 
that contain the same self-energy piece as diagram 2(a) also contain the In Q-contribution 
that is characteristic for this diagrammatic element. The results are again given in Table |. 

Table | lists the contributions of all diagrams to the coefficient a 2 in Eq. (|2.5|) . In 
some cases we have found it convenient to combine some diagrams before evaluation of the 
integrals, and this is indicated in the table. The coefficients a\ and a 2 i og have been calculated 
in Ref. [TBI with the result, 



a x = -4tt/3 



09 



(vr 2 - 4)/2 



Table | contains three integrals which we could not reduce to tabulated ones, viz. 



h 



1 dx 
x 



In 



x 



1 + x 



4.207. 



(2.12) 



(2.13a) 



dx x 



In 



1 — x 

1 + x 



2.772 . 



(2.13b) 



4 r°° dx . 
1 3 — — / — (arctanx 

71 JO X 2 



1 

1 arctan x 

x 



7.716.. 



(2.13c) 



Summing all the contributions listed in Table [I] we obtain 
55 

a 2 = A\n(Q/k F ) + — vr 2 - 14 In 2 + 7 - h + I 2 - h = Aln{Q / k F ) + 3.22... . (2.14) 
36 

This concludes our calculation of the zero-temperature conductivity. In order to compare 
with, and make explicit predictions for, experiments it is desirable to convert this result into 
the corresponding one for the mobility at nonzero temperature. In order to do this, we will 
need the density of states to second order in the impurity density as well as the conductivity. 
The density of states is easily obtained from the Green's function via the relation, 



N(e) = — £!m (<#> = 0) 



7T 



k,p 



(2.15) 



and to second order it is sufficient to consider the diagrams shown in Fig. The calculation 
is easy, and we obtain, 



N(e) = iV (0) 



1 + 



7 



2 ]n(Q/2y/e) -2 In 2) + 0( 7 



,2eJ V 2 

Here = 2N F = mkp/7c 2 = m^fejix 2 is the free electron density of states. 



(2.16) 
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C. The Electron Mobility 



In a time-of-flight experiment like the one described in Ref . O the measured observable is 
the electron mobility, which is given by /i(T) = a(T)/en(T), where cr(T) is the temperature 
dependent conductivity, and n(T) is the electron particle number density. cr(T) can be 
obtained from Eq. (|2.5| ) by means of the Kubo-Greenwood formula,!! and n(T) from the 
density of states, Eq. (|2.16|) . We thus have, 



de 



-df 

de 



cr e 



def(e)N(e) 



(2.17) 



with o-(e) and N(e) from Eqs. (p3|), ( 01% ), ( gig) , and (|2~T6| ), respectively, and /(e) the 
Fermi function. In writing all of these quantities as functions of e = k 2 F one has to keep in 
mind that the electronic mean free path / is the energy independent parameter which has 
to be kept fixed. We are interested in the limit of small electron density, where the Fermi 
function can be replaced by a Boltzmann distribution, /(e) ~ exp (/?//) exp(—(3e). Doing the 
integrals we obtain Eq. ( |1.3a| ) for the temperature dependent mobility with Hi and \i2iog as 
given by Eqs. fll.3bUl.3q 1 ), and 



A*2 



1 

16 



7l~ 

36 



(55 + 9C) - C + 8 - 10 In 2 - h + h - h 



-2/z 2W ln2 = 0.236... 



(2.18) 



where C is Euler's constant. As expected (see the discussion after Eq. ( [2.11| )) the mobility 
is independent of the cutoff Q. 



III. DISCUSSION 

Since the mobility is directly measured in a time-of-flight experiment of the type reported 
in Ref. [T3|, the density expansion for the mobility as given by Eqs. ( |1.3j ), ( |2.1£j| ) can be 
directly compared with experiment. We split the discussion of the experimental relevance 
of our result into three separate questions: (1) How accurately does our quantum Lorentz 
model describe electrons injected into Helium gas? (2) How do our results compare with 
existing data? (3) What kind of experimental effort would be necessary to unambiguously 
detect the logarithmic term in the density expansion? 

The first question has to be divided into idealizations that are inherent in the model, 
and effects of additional approximations we have made. Let us start out with the latter. 
As mentioned in the Introduction, there are three independent length scales in the model, 
viz. the mean scatterer separation n _1//3 , the mean-free path /, and the thermal wavelength 
A. Consequently, one can form four different dimensionless densities, viz. nX 3 , n\ 2 a s , nXa 2 , 
and na 3 . The first two do not appear in the Lorentz model.0 Of the other two, the first one 
is essentially the ratio of the thermal wavelength to the mean-free path, while the second 
one describes an excluded volume effect that is also present in classical systems. We have 
kept only the leading (for X/a s » 1) one of these two parameters, nXa 2 s . That is, we 
have neglected terms of relative order a s j X. At Helium temperature, and with the electron- 
Helium scattering length [a s = 0.63A), the value of this small parameter is a s /X ~ 10~ 3 . 
The prefactor can be estimated from classical Enskog theory, or from the exactly known 
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result for the classical Lorentz model,!!, and turns out to be of order unity. On the other 
hand, a typical value for the expansion parameter x m Eq. (ECU) is X ~ 0.1.0 In our 
expansion, excluded volume corrections to the term linear in x w iU therefore roughly be of 
the same order as terms ~ x 3 - We can thus safely neglect excluded volume effects, unless 
X is chosen to be too small, see below. It is also worthwhile to point out that the excluded 
volume terms have a different temperature dependence than the leading terms that were 
kept in Eq. ( |1.3a| ). In principle it would therefore be possible to separate these contributions 
experimentally, although this may be hard to do in practice. 

Another approximation has been our restriction to pure s-wave scattering. This amounts 
to neglecting corrections of 0{kp/Q) ~ 0{k F a s ) in Eq. (|275|) , which translates into correc- 
tions of 0(a s /X) in Eq. flOaD . Non- s-wave scattering effects and excluded volume effects 
are therefore comparable. 

Idealizations that are inherent in the model include the static nature of the scatterers, the 
assumption that the scatterers are uncorrelated, and the single-electron approximation. The 
effects of static correlations between the scatterers can be easily estimated from, e.g., Baym's 
formula for the electron scattering rate, which can be obtained as a variational solution of 
the Boltzmann equation.^ For the case of static impurities, Baym's formula contains a static 
impurity structure factor S(q) under the momentum integral that determines the scattering 
rate. For uncorrelated scatterers, S(q) = 1. Deviations of S(q) from unity, and hence 
static correlations between the scatterers, are again due to excluded volume effects which 
we have estimated above. The dynamics of the Helium gas lead, at temperatures around 
4K, to corrections of the same order, which can be seen as follows. A typical value for the 
scattering time is r « 10~ 12 s. During this time the thermal velocity of the Helium atoms 
leads to a displacement d on the order of 1A. This is comparable with the scattering length, 
and the effect on the electron mobility is thus on the order of d/X ~ a s /X, which is again 
the magnitude of the excluded volume effects. Finally, one has to estimate the effects of the 
Coulomb interaction, since any real experiment has to deal with more than one electron. Let 
us use the parameters of Schwarz's experiment^ for this purpose. A typical value for the 
electron current density was j ~ 10 -12 Ccm~ 2 s~ l , and for the drift velocity v ~ 10 4 cm/s. 
This corresponds to an electron density n e m 10 3 cm~ 3 . The Coulomb energy E c m e 2 /n~ 1//3 
is then less than one percent of the kinetic energy, ksT. If desirable, Coulomb effects can 
be made even smaller by decreasing the electron density. 

We conclude from the preceding discussion that both the model and our additional 
approximations are adequate for a description of electrons in Helium gas in a parameter 
range of interest to us, namely at temperatures and densities such that our expansion, 
Eq. ( |1.3a| ), is meaningful. Let us now discuss the relation of our result to existingEl and 
possible future experiments. The main remaining uncertainties in this relation arise from 
the terms of o(x 2 ) and 0(xa s /X) in Eq. ( |1.3a| ). The former are undoubtedly nonanalytic, 
but neither their functional form nor their magnitude are known.ll In order to estimate their 
importance we therefore have to rely on some assumptions. First, it is plausible to neglect 
the nonanalyticity, and to assume that the coefficients in the x-expansion are all of roughly 
the same magnitude. In the Kubo-Greenwood integration, Eq. (|2.17f) , the x 3 -term picks 
up an extra factor of t/tt, and so we estimate o(% 2 ) ~ fi 3 x 3 with — 2 A /7r / u 2 ^ 7/3 ^ 2^fnji 2 - 
Here we have neglected the, presumably weak, x-dependence of /x 3 , and have allowed for 
a safety margin in the form of an extra factor of 2. Obviously, the relative effects of [13, 
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and the other unknown terms of higher order, become smaller with decreasing values of \- 
On the other hand, the excluded volume effects are of 0(xo s /A), and become important if 
X becomes very small. However, as long as they are small compared to the second order 
terms which we keep in Eq. (|1.3a| ), i.e. as long as \ 2 > X a s/^, or \ > a s/^, they can be 
neglected. This means there is a window of ^-values for which the excluded volume terms 
are negligible, but the higher order in x-terms are not yet important. 

We now define, as a convenient quantity directly comparable with experiment, 

fix) = W)/Hb - 1 - mx] /X 2 ■ (3-la) 
Our theoretical prediction for this quantity is, 

fix) = V-ilog In X + V-i ± /^22\Arx , (3.1b) 

with /i 2 from Eq. ( |2.18|) , and \in og from Eq. (|1.3cj ). Let us consider the experimental results 
obtained by SchwarzB at Helium temperature. At T = 4.2K, a He gas density n = 10 21 cm~ 3 
corresponds to x = 1, and data were obtained for x as l°w as 0.08. In Fig. [7] we show the 
theoretical prediction, Eqs. ( [3.1|) , for < x < 0.7 together with Schwarz's data. The error 
bars shown assume a total error of 3% in \ij \ib an d 4% in \. To illustrate the effect of the 
logarithmic term the figure also shows what the theoretical prediction would be if Oiiog i n 
Eq. ( |2.5| ) was zero. 

We are now in a position to draw the following conclusions. (1) The existing data are 
certainly consistent with the existence of the logarithmic term, but are not accurate enough 
to be conclusive. (2) A repetition of this experiment in the region 0.1 < x < 0.2 with 
an accuracy improved by at least a factor of 10 would be sufficient for a convincing test 
of the logarithmic term's existence. This range of ^-values is particularly suitable because 
excluded volume effects are negligible. At larger %-values the uncertainty due to the x 3 -terms 
makes the theoretical prediction meaningless, and at lower %-values errors in determining x 
translate into very large errors in f(x)- Also, the excluded volume effects become noticable 
at smaller x- (3) The type of experiment discussed, i.e. a time-of-flight measurement for 
electrons injected into Helium gas, probably constitutes the most favorable opportunity for 
an experimental check of the existence or otherwise of the logarithmic term in the density 
expansion for transport coefficients. We hope that the foregoing discussion will stimulate a 
new precision experiment on this system. 
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APPENDIX A: EVALUATION OF THE INTEGRALS J ++ (q) AND J+ (q) 

We have seen in Section II that the diagrams can be conveniently expressed in terms 
of the two integrals J ++ (q) and J + ~{q) given by Eqs. (|2.10|) . In Refs. p~5|JT6| approximate 
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representations for these integrals were given, which are not sufficient for the purposes of 
the present paper. We therefore evaluate J ++ (q) and J + ~(q) exactly. The method we are 
using actually allows us to calculate two slightly more complicated functions, viz. 



and 



J ++ (g,c)= J 
J + -(q,uj) = J 



dk 



dk 



e — k 2 + ry e + ui — (k — q) 2 + ry 



e — k 2 + 27 e + u — (k — q) 2 — ry 



(Ala) 



(Alb) 



J ++ {q) and J + ~{q) are obtained as the values of these integrals at oo = 0. 

To proceed we introduce spherical coordinates (k, 9, 0). The integral over the azimuthal 
angle <fi is trivial and gives a factor of 2-7T. The integration over the polar angle 9 e [0, ir] 
can be written as an integral over x = cos 9 extending from -1 to +1. Using the symmetry 
of the integrand then allows us to write 



J +U (q,oo) = 2tt f 1 dx f 
Jo Jo 



dk 



k 2 



a 



k 2 



1 



+ 



1 



V - k 2 - 2kqx f3 v -k 2 + 2kqx 



(A2) 



where a = e + ry, j3 v = e + oo + ivy — q 2 and v — ±1. Since the integrand is an even function 
of k, one can now extend the /^-integration to the interval [— oo, oo], and evaluate the integral 
by means of the residue theorem. The remaining integrals over x are elementary, and one 
finds, 



27T 



In 



In 



'2vy/aj3 v + q(a - f3 u 



a 



Vq 2 + Pu + q 



2v^t(5 v - q(a - fa) + (a + Pv)\/q 2 + f3 V/ 
After some simplifications we obtain, 

r2 



(A3) 



J ++ (q,oo) = 



17T 



In 



y/e + oo + ry + y/e + iy - q 



(A4a) 



and, 



ITT 



J + (q, oo) = In 



\Je + oo - iy - \[e + ry - g 



_ Ve + ^ - *7 - V e + ^7 + <? / 
In the limit u; — > 0, J + ~ can be further simplified, and we finally obtain, 



(A4b) 



and 



J ++ (q) 



In 



'2y/e + vy + q* 
2y/e+lrf-q / 



2tt 2 



arctan 



2lmy/e + iy 



(A5a) 



(A5b) 
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For small 7 these expressions agree with those obtained previously.lll 

For explicit perturbative calculations the following formal small-7 expansions are useful, 



J++(g) 



3 2 

— Q(q-2yft) -1— In 
q q 



q-2^l 



q + 2^Te 



e I 4e - q 2 



i n 3 y^5(q-2V~e))+0( 1 2 ) 



7i 

Q 



7 2vr 2 
v/e g 2 



0( 7 2 ) 



(A6a) 
(A6b) 



APPENDIX B: EVALUATION OF DIAGRAM 2(a) 

To calculate the contribution of diagram 2(a) to the zero-temperature conductivity we 
use the expressions for J ++ (g) and J + ~(q) from Appendix [S]in Eq. (j2.8|), which we rewrite 

as, 

a W = ° B {£) ^72" [ Re ( " Z7) jC*^ ( J++ ^)) 2 

+Re [ Q dq q 2 J ++ (q) J + ~(q) - e Re / Q dg g 2 J ++ (g) J + ~(<?) 

Jo JO 



-e 7 Re [ Q dqq 2 J ++ (q)-^-J ++ (q) + eRe [ Q dq q 2 fj ++ (g)" 2 
Jo cry Jo v 



+ 0( 7 2 ) • (Bl) 



o?7 

Because there is a factor 7/2e multiplying the integrals it is sufficient to evaluate the latter 
to linear order in 7. Of the five integrals the first two can be straightforwardly evaluated 
by using the expansion given in Eq. ( |A6| ) . The same is true for the third integral, although 
the expansion procedure leads to an ill-defined integral ^dx/{x — 1), which has to be 
interpreted in a principal values sense (so it is zero). Alternatively, one can use the exact 
expressions for J ++ and J + ~ from Eqs. ( |A5| ) , write the arctan as an auxiliary integral, 
arctanx = x Jq dy/(l + y 2 x 2 ) } and use complex analysis. The result is the same, viz. 

Re J® dq q 2 J ++ (q) J + -(q) = tt 6 (Q - 2Re^/e + 27) - (7 4^/^) ln(Q/2^) + 0( 7 2 ) , 

(B2) 

where we have dropped terms that vanish for Q — > 00. The fourth term requires the 
evaluation of an integral, 

J = lim Im [ dq q In ( —, '— - ) , , (B3a) 

s->o J-Q y2^/e + vy -qj q + 2y/TTry + 5 q - 2 y /TT 7 H ~ & 

where we have made use of the symmetry of the integrand, and have shifted the poles off 
the branch cuts of the logarithm. Standard complex analysis techniques yield, in the limit 
of large Q, 

J = -2tt [ln(Q/2Ve) - In 2] + 0( 7 2 ) . (B3b) 

Finally, the fifth integral can be related to the fourth one by Taylor expanding in 7. Adding 
all contributions one obtains entry 2(a) in Table |. 
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FIGURES 



FIG. 1. The simple bubble (a), and two diagrams containing the 'triangulated' Green's function 
(b), (c). Together with diagrams (b) and (c), their complex conjugates (c.c.) also contribute. 



FIG. 2. The diagrams that were considered previously in Refs. 15,16. 
FIG. 3. Infinite resummations that are derivatives of diagrams 2(d)-2(h). 
FIG. 4. Skeleton diagrams with three impurity lines 
FIG. 5. Skeleton diagrams with four impurity lines 
FIG. 6. The Green's function to second order in the impurity density 



FIG. 7. The reduced mobility /, as defined in Eq. ( 3.1a| ), vs. the density parameter x = A/7rZ. 



The theoretical prediction is for / to lie between the two solid lines. The experimental data are 
from Fig. 9 of Ref. |l4| with error bars estimated as described in the text. The broken lines show 
what the theoretical prediction would be in the absence of the logarithmic term in the density 
expansion. 
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TABLES 



TABLE I. Values of the diagrams as shown in Figs. Only the contribution to the coefficient 
o"2 in Eq. (2.5) is given, for contributions to lower orders see Refs. 15,16. See the text for further 
explanation 



diagram 


C2 


diagram 


0"2 


1 (a) 


3/2 


4 (f) 


-h + 2/ 2 


1 (b) 


-7/2 


4 (g) 


41n2 - 2 


1 (c) 


1 


4 (h) 


2-41n2 


2 (a) 


61nQ/^ + 4- 12 In 2 


4 (i) 


A/2 - h 


2 (b) 


8 


5 (a)+(b) 





2 (c) 


-2-41n2 


5 (c)+(d) 





2 (d) 


7T 2 /2 + 3/l/2 


5(e) 


7T 2 /9 


2 ( e )+(f)+(g)+(h) 


-1 -vr 2 - 2 In 2 


5(f) 





3 (a)+(b)+(c) 


vr 2 /4-21n2 


5(g) 


2vr 2 /3 


3 (d) 


-h 


5 (h)+(i) 


vr 2 


4 (a) 


-ii/2 


5 (j)+(k) 





4 (b) 


-Ii/2 


5 (l) + (m) 





4 (c) 


-/1 


5 (n) + (o) 





4 (d) 


-21n(Q/Ve) +4 In 2 


5(P) 





4 (e) 


2 In 2 - 1 
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k k-q k 




1 



